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The role of spherical quantum shells in the competition between fusion and quasi-fission is studied for reactions forming heavy 
■elements. Measurements of fission fragment mass distributions for different reactions leading to similar compound nuclei have 
'been made near the fusion barrier In general, more quasi-fission is observed for reactions with non-magic nuclei. However, 
■the '^^Ca-i-^^^Pb reaction is an exception, showing strong evidence for quasi-fission, though both nuclei are doubly magic. Time- 
[dependent Hartree-Fock calculations predict fast equilibration of NjZ in the two fragments early in the collision. This transfer 
"of nucleons breaks the shell effect, causing this reaction to behave more like a non-magic one in the competition between fusion 
^and quasi-fission. Future measurements of fission in reactions with exotic beams should be able to test this idea with larger NIZ 
■asymmetries. 
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, Quantum shell effects play a key role in the structure and 
"stability of atomic nuclei, as they do in the periodic chemical 
.properties of the elements. Where there is a large energy gap to 
■the next quantum level, the total number of protons or neutrons 
.filling all levels below the gap is referred to as a magic number 
Tn particular, magic nuclei have a smaller mass per nucleon than 
.their neighbours. The variation of the magic numbers across the 
■nuclear chart is crucial to build our understanding of the nuclear 
.quantum many-body system. One major challenge is to define 
■the magic numbers in the region of the superheavy elements 
^(SHE), with Z > 110 protons H H S H. Associated with 
this, atom-by-atom measurements of the chemical properties of 
.SHE are testing the predicted strong relativistic effects on the 
■electrons which modify the periodic Table \3\ . 
. SHE up to Z - 118 have been synthesised in fusion reac- 
■fions of heavy nuclei, either using ^^^Pb and ^°^Bi targets 
,or Ca beams on actinide targets Production cross sec- 

tions are, however, extremely small (of the order of a few pico- 
.barns), and a good understanding of the reaction mechanisms 
■is needed to optimise their production. To achieve a compre- 
hensive global picture of SHE formation is very challenging, as 
many variables may affect fusion probabilities. These include 
collision energy, mass-asymmetry, deformation and orientation, 
isospin, and shell structure of the colliding nuclei. These vari- 
ables are often strongly entangled, making it difficult to isolate 
the effect of a single variable. Furthermore, these properties 
evolve dynamically, thus it is necessary to understand the dif- 
ferent associated time scales. 

The early stage of the collision is a crucial step in SHE for- 
mation, where the initial conditions are the most important. 



These determine the configuration at which the colliding part- 
ners have dissipated their kinetic energy, thus determining the 
shape of the initial fragile dinucleus. This can break apart, gen- 
erally after multiple nucleon transfers (mainly from the heavy to 
the light partner), in a process called quasi-fission (QF) iSlllS]- 
Alternatively it may reach compact shapes, fusing to form a hot 
compound nucleus (CN), which can lead to formation of a SHE 
through neutron evaporation in competition with CN fission. 
Although the CN survival probability against fission is very 
small, its decay width is governed by the well-known equations 
for statistical decay, which should allow prediction of the rela- 
tive survival probabilities from different fusion reactions. This 
is not the case for QF, which is a completely dynamical pro- 
cess, and depends on many variables which can be different for 
different reactions. The nature of the two fission processes are 
reflected in their time scales, which can be very different. Typ- 
ical times scales for QF are shorter than 10"^" s [tI HImIj 
but can be longer than 10 s ifioll for fusion-fission. To ef- 
ficiently form SHE, the entrance channel conditions should be 
chosen to minimise the QF probability, which is dominant in 
reactions forming SHE. Beyond the basic principle of minimis- 
ing the Coulomb energy in the entrance channel, a quantitative 
understanding of the effects of nuclear structure on the compe- 
tition between fusion and QF is a key missing ingredient. 

At collision energies above the fusion barrier, a systematic 
analysis showed that closed shells in the colliding nuclei have 
a relatively small effect on fusion probabilities 1 1 1]. However, 
at energies around the barrier, the competition between fusion 
and QF is known to be affected by (shell-driven) nuclear defor- 
mation and orientation IH, [H [3, [Bl E [13, [Hi ■ Spherical 
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shells may also be important, resulting in so-called "cold val- 
leys" in the potential e nergy surface, which lead to the compact 



CN configuration 
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Figure 1: (colour online). Measured mass-angle distributions for each reaction (upper panels). The factor multiplies the maximum counts of the logarithmic colour 
scale (right). In the projected mass ratio spectra (lower panels) the scale factor multiplies the counts scale on the left. The difference between the scale factors is due 
to the various statistics obtained with each system. Gaussian fits to the region around Mn=Q.5 are shown (turquoise lines), whose standard deviations cr are given in 
Tabled Gaussian functions with ctur = 0.07 (thin red lines) are shown for reference. 



leys may also be favoured because energy dissipation should be 
weaker, allowing greater inter-penetration before the initial ki- 
netic energy is dissipated 23']. These effects may be vital 
in the recent successful synthesis of SHE |1, 2, 3,^. However, 
the interplay of spherical shells with other degrees of freedom, 
such as the isospin of the two colliding nuclei, has not yet been 
investigated. 

In this letter, the role of spherical shells (magicity) on the QF 
probability is first demonstrated through fission measurements 
for reactions with relatively small initial isospin asymmetry, or 
more precisely N/Z asymmetry, quantified by the difference be- 
tween the N/Z ratios of the initial colliding nuclei A(A^/Z),. 
Then, the case of a magic reaction with large A(A^/Z), is in- 
vestigated. The time scales for QF and isospin equilibration are 
investigated with the help of calculations, and used to explain 
the measurements in terms of the dynamical interplay between 
isospin asymmetry and spherical shells. 

Measurements were made using the 14UD electrostatic ac- 
celerator at the Australian National University. Pulsed beams 
of 111 MeV i^O and 213.5 MeV ^^^Ca, and DC beams (giv- 
ing higher intensities) of 212 MeV '^Cei, 213 MeV "^^^Ca, and 
230, 235 MeV '^'^Ti were produced from metallic ""'Ca and ™'Ti 
samples. Isotopically enriched targets of ^"'^Pb (420 fig/cm^ 
self-supporting), and ^o^PbS (30 yug/cm^), ^ooHg (15 fig/cm^) 
and ^^^UF4 (400 /vg/cm^), evaporated onto ~15 fig/cm^ ™'C 
backings, were mounted on a target ladder whose normal was 
at 60 ° to the beam. Binary reaction products were detected 
in coincidence using two 28x36 cm^ position-sensitive mul- 
tiwire proportional counters on opposite sides of the beam, 
covering laboratory scattering angles of 5 ° < < 80 ° and 
50 ° < 9 < 125 °. For the pulsed beams, the measured positions 
and times-of-flight allowed direct reconstruction of the frag- 
ment velocities il4ll . With DC beams, the velocities were deter- 
mined from the time difference between the two fragments fl^. 



assuming binary reactions and full momentum transfer, which 
will be valid for the low fissility targets used il4ll . Following 
iterative correction for energy loss in the target, the mass ratio 
Mr - mi /(nil + ni2) (where nii and m2 are the two fragment 
masses) and the centre-of-mass (cm.) scattering angle Oc.m. 
were deduced. Since both fragments are detected, the mass- 
angle distribution (MAD) is populated twice 1,24,1 . at {Mr, Oc.m) 
and(l - MR,n - Ocm). 

The MAD for the reactions measured are shown in the up- 
per panels of Fig. [1] The reactions with Ca and Ti beams form 
isotopes of the elements No (Z=102) and Rf (Z=104), and in- 
volve similar charge products in the entrance-channel. The 
i6q^238u reaction forms Fm (Z=100), but with less than half 
the entrance-channel charge product. In the measurement the 
azimuthal coincidence coverage was essentially 90 ° for all 6, 
thus the number of events in each MAD bin is proportional 
to the angular differential cross section dcr/dflf ,„ . Note, how- 
ever, that every MAD has a different coefficient of proportion- 
ality due to the varying statistics obtained for each reaction. 
The intense bands at extreme Mr values correspond to elas- 
tic and quasi-elastic (QE) scattering, while fission-like events, 
associated with either fusion-fission or QF, are spread around 
Mr - 0.5. Note that, in our measurements with Ca and Ti 
beams, both fusion and QF occur at similar partial waves. In- 
deed, the beam energies correspond to below-barrier energies, 
as can be seen from Table[T] where centre-of-mass energies and 
calculated barrier energies are given. Thus the angular mo- 
menta involved are low, and those of fusion and quasi-fission 
are bound to show a large overlap. 

For the heavier projectiles, the fission-like events clearly 
show a correlation of fragment mass with angle, resulting from 
the short reaction times (< 10"^" s) [2, 8, 9]. For example, for 
the '*'*Ca-i-2°'*Pb reaction, the Mr centroid for 125 ° < Ocm. < 
135 ° is 0.511 + 0.004. Although the deviation from Mr=0.5 
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is small, it is much larger than the statistical uncertainty, and 
as clearly seen in the MAD, varies consistently with Oc.m.- Ref- 
erence measurements for the reactions of "'O with ^*'**Pb, and 
with ^^^U (shown in the left-most panel of Fig. [1) give essen- 
tially no correlation of mass with angle, consistent with much 
longer fission times. 

The lower panels of Fig. [T] show the M« projections of the 
MAD spectra above. The widths of the fission-like fragment 
mass distributions are expected to be larger in the presence of 
QF than if only fusion-fission is present [S] 17, 250 ■ To char- 
acterise the Mr distributions for the fission-like events, and 
to allow comparison with previous work \Tf\, they were fit- 
ted with Gaussian functions, within the range 0.34<Mr<0.66 
(turquoise curves in Fig.[TJ so as to exclude deep-inelastic and 
QE events. For ^''O+^'^^X], we choose 0.2<Mr<0.8 as only 
fission-like events were detected. The fitted standard deviations 
a MR are given in Table [T] together with the value for 218 MeV 
'*^Ca-i-2''^Pb from Ref. [27]. Since it may well be that the true 
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distributions are not single Gaussians 
we also compute the standard deviation I^mr of the data points 
in the same 0.34<Ms<0.66 range, which are also given in Ta- 
ble[T] Of course, <jmr and "Lmr are different quantities with dif- 
ferent values, but they both constitute a measure of the width of 
the fission-like fragment mass distributions, the latter indepen- 
dent of any assumed shape. As will be seen, the two quantities 
do exhibit the same trends, and together with the reasonable 
reproduction of the experimental data by the Gaussian fits, sug- 
gest that the fitted standard deviations ctmr give a reasonable 
representation of the mass width of the fission-Uke events, with 
a single parameter. 

In order to investigate the influence on quasi-fission of spher- 
ical shells in the entrance channel, we plot in Fig. [2] the widths 
{(Tmr and 1.mr) of the fission-like fragment distributions as a 
function of the number A^,„ of entrance channel magic numbers 
(given in Table[T]). The possible proton and neutron magic num- 
bers for projectile and target nuclei are Z - 20 and = 20, 28, 
and Z = 82 and = 126 respectively. An upper limit to 
the standard deviation ctmr for fusion-fission (o-fiss) of 0.07- 
0.08 can be taken from the present and previous |@, [l3l mea- 
surements for '^O-H^^^U. This is only an upper limit as it was 
shown lfl4ll that QF contributes to fission-like events even in 
this reaction. This range is indicated by the horizontal band in 
Fig. ma). Only the ^^Ca-i-^o^.zosp^ ^j^^^ jjg jj^jg ^^^^^ ^jj 

other reactions have larger widths, indicating the presence of 
QF |H,[l3>|25t], a result consistent with the observation of a de- 
pendence of mean fragment mass with angle in the measured 
MAD presented in Fig.[T] 

Apart from the ^''CaH-^''^Pb reaction, discussed later, a clear 
correlation is seen in Fig. |2] between the entrance channel 
magicity, quantified by N„,, and the amount of QF, related to 
cmr - o-fiss- the less entrance-channel magicity, the more QF. 
As discussed in the introduction, this correlation could result 
from cold valleys in the potential energy surface CM Edi El 



2211 and/or weaker energy dissipation 01 1112311 . both eff'ects be- 
ing associated with the spherical shells. As a result, a greater 
inter-penetration of the two nuclei should then be achieved, 
leading to a higher fusion probability, and, consequently, a 
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Figure 2: (colour online), (a) Standai'd deviations ctmr of the Gaussian fits 
of the fission-like fragment mass distributions as a function of the number of 
magic numbers in the entrance channel N,„. The horizontal band shows the 
upper limit of o"mk for pure fusion-fission (i.e., without QF). (b) Standard de- 
viation l^MR in the 0.34<Mr<0.66 range. The dotted line coiresponds to a flat 
distribution. When not shown, the statistical uncertainties are smaller than the 
size of the points. The dashed lines guide the eye. 
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Table 1: Standard deviation 'Zmr of fission-like fragment mass distributions and standard deviation (Tmr of their Gaussian fits (see text) for each reaction, with 
statistical uncertainties. Ni„ is the total number of inagic numbers in target and projectile, and A(A'/Z), is the difi"erence between their ^ ratios. Centre-of-mass 
energies i?r and theoretical barriers B,i, from the proximity model | 30] are in MeV. 
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254pjj^ 
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80.3 
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0.59 


0.081 +0.001 


0.073 + 0.001 


Present work 


48Ti+200Hg 


248^0 


185.5 


190.9 





0.32 


0.237 + 0.025 


0.090 + 0.001 


Present work 


4STi+2"8pb 


256p,f 


190.9 


194.4 
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0.35 


0.121+0.004 


0.082 + 0.001 


Present work 


44Ca+204pb 


248^0 


174.4 


178.0 
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0.29 


0.114 + 0.002 


0.081 +0.001 


Present work 


48ca+204pb 


252No 


172.4 


175.8 
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0.09 


0.084 + 0.008 


0.073 + 0.004 


Present work 


40ca+208pb 


248^0 


179.1 


179.5 
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0.54 


0.126 + 0.004 


0.083 + 0.001 


Present work 


48ca+208pb 


256No 


177.1 


175.0 


4 


0.14 


0.068 + 0.002 


0.064 + 0.002 


From Ref. [27] 



smaller QF probability. This interpretation is also supported 
by the observation of relatively high fusion-evaporation cross- 
sections (up to ~ 3 i-ib) in the "^^^Ca+^^^^Pb system as compared 
to reactions with non-magic targets with similar masses \2^. 

Before accepting this conclusion, the possible effects of a 
number of additional variables should be considered. The com- 
parison of the widths is strictly valid for reactions forming 
the same CN (here, the three reactions forming ^'^'^No), under 
the same conditions. The known dependence of the standard 
deviation ctmr on excitation energy for these reactions is too 
weak lfl4ll27ll . and the difference in energies too small (e.g., ex- 
citation energies in ^^Ti+^^^Hg and ^'^Ca+^^^Pb differ by only 
0.4 MeV) for differences in excitation energy to affect the con- 
clusions. The '*'^Ti+^*'**Pb reaction has the largest entrance- 
channel charge product, and forms the heaviest and most fissile 
nucleus, thus without shell effects, the largest standard devia- 
tion (Tmk might be expected. This is not what is observed, so 
we conclude that the large changes in ctmr must be related to 
the differing magicity in the entrance channel, rather than prop- 
erties of the composite system. 

There is one reaction that does not follow the systematic be- 
haviour shown by the others, namely '"^Ca+^^^^Pb. Fig.|2]shows 
clearly that it demonstrates strong evidence for QF (ctmr - 
0.13), despite having maximal magicity Nm - 4. We propose an 
explanation below which does not invalidate the link between 
magicity and QF probability seen for the other reactions. To 
solve this puzzle, it is sufficient to invoke the fast isospin equili- 
bration resulting from nucleon transfer Detailed measurements 
of reaction product yields 113 ill and angular distributions |32] 
have shown that systems with strong isospin asymmetry in the 
entrance channel (like '**'Ca+^°**Pb |32]) undergo a rapid (al- 
though incomplete) isospin equilibration in the early stage of 
the collision, through the transfer of nucleons ll32ll . 

The time-dependent Hartree-Fock (TDHF) theory has suc- 
cessfully described transfer in N/Z asymmetric reactions (for 
example Refs. ||33, 34, 35l[3^). Here, it is used to investigate 
the timescale of isospin equilibration via transfer. In TDHF, 
each particle evolves independently in the mean-field generated 
by all the others. The TDHF formalism is optimised for the pre- 
diction of expectation values of one-body operators, such as the 
average N/Z ratio in the fragments. The tdhf3d code is used 



with the SLy4£/ parameterisation of the Skyrme functional 137 1 . 
The TDHF equation is solved iteratively in time, with a time 



step Af = 1.5 X 10"-"* s, on a spatial grid of 56 x 56 x 28/2 
points with a plane of symmetry (the collision plane), and a 
mesh size Ax = 0.8 fm (see |38] for more details). The initial 
distance between the nuclei is 22.4 fm. 

The results of the TDHF calculations of N/Z equilibration 
between two colliding nuclei are shown in Fig. |3] as a func- 
tion of their contact time. The difference in the N/Z ratios of 
the two nuclei before any transfer of nucleons is denoted by 
A(N/Z)j, shown by the full circles in Fig. [3] and also given in 
Table [1] The curves show the calculated evolution of the differ- 
ence between the N/Z ratios of the outgoing (final) fragments 
(A(N/Z) f) for each reaction, as a function of contact time, de- 
fined as the time during which the neck density exceeds half 
the saturation density po/2 - 0.08 fm"-'. The contact time is 
varied by making calculations at angular momentum Lh from 
~ 20h to IQh. The energies of the collisions are the same as 
in the experiment. For L < 20, most of the systems undergo 
capture resulting in fusion, whose timescales are too long for 
the TDHF calculations, or strongly damped collisions. For 
the smallest contact times (associated with large L), the nu- 
clei scatter (in)elastically and no change in isospin occurs (i.e., 
A(N/Z)j ^ A{N/Z)j) as seen in Fig. |3] For the 48Ca+2"4.208pb 
reactions, the initial isospin asymmetry is small, and no change 
in isospin occurs with increasing contact times. The fact that, 
for these reactions, A(N/Z) j never reaches zero is typical for 
mass asymmetric reactions [39]. For the other reactions, as 
the contact time increases, the A(N/Z) j approaches the same 
isospin asymmetry. In particular the most N/Z asymmetric re- 
action, ''"Ca+^°**Pb undergoes a large reduction of A(N/Z)f, in 
agreement with experiment 132]. Using a particle number pro- 
jection technique ]40], the most probable outcome for this re- 
action after a contact time of ~ 2.7 x 10"^' s (calculated for 
L = 20), is found to be "^^Ar+^^^Po. This calculation also 
gives the probability of remaining in the entrance channel (and 
thus of conserving its entrance channel magicity), which is 
Paz=oP&n=o ^ 0.083 X 0.002 ^ 1.7 x 10""*, a negligible proba- 
bility. However, for the '***Ca+^''**Pb reaction, even for a contact 
time as long as ~ 3.5 x 10"^' s (not shown in Fig|3]l, this prob- 
ability is still 0.76 X 0.57 ^ 0.43, giving a much larger survival 
probability for the initial magic numbers in this reaction. 

If the nucleons are transferred and the magicity is lost early in 
the collision, the system should behave more like a non magic 
system. On the contrary, if isospin equilibration takes place 
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Figure 3: (colour online). Final NjZ asymmetry of the fragments (lines) as 
a function of contact time of the fragments in zeptoseconds (1 zs=10""' s), 
defined as the time during which the neck density exceeds half the saturation 
density poll = 0.08 fm"^^. The initial values A(N/Z)i are shown by full circles. 
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Figure 4: (colour online). Simulated mass-angle distributions (middle panels) 
for '"'Ca+-°*Pb following |9], for the quasi-fission time distributions shown in 
the upper panels. The average time (t) of these distributions is indicated. The 
projections on the mass ratio axis are shown in the lower panels. 



on a time scale similar to that of QF, then the magicity in the 
entrance channel could still significantly enhance fusion. Ac- 
cording to the TDHF calculations (Fig.[3]l, ^°Ca-i-^°'*Pb experi- 
ences a high degree of isospin equilibration for sticking times 
< 2 X 10"^' s. This is in agreement with experimental observa- 
tions II32II of a high degree of N/Z equilibration in deep-inelastic 
collisions before many nucleons have been exchanged. This 
time has to be compared with the typical time scale for QF. Mi- 
croscopic quantum theories cannot yet model such collisions 
from first principles 1,41.1 . Thus, to obtain the QF time for the 
reactions studied, MAD have been simulated using a classical 
trajectory model [9]. MAD were calculated for three differ- 
ent QF time distributions shown in the upper panels in Fig. H] 
whose mean times varied from 3.5x10"^' s to 14x10"^' s. The 
calculated MADs corresponding to these mean times are shown 
in the middle panels of Fig. |4l whilst the bottom panels show 
the predicted mass ratio spectra. The shape of the experimen- 
tal data (Fig.[T]right panel) is best reproduced with an average 
time scale of 14x10"^' s, which is much longer than the time 
for isospin equiUbration. Isospin equilibration leading to loss of 
magicity occurs early in the ""'Ca+^^^^Pb collision, which thus 
may be expected to exhibit QF properties closer to non magic 
systems. This is what is seen experimentally, as clearly shown 
inFig.|2] 

Finally, let us note that previous measurements 
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of excitation functions for capture reactions (including both 
fusion-fission and quasi-fission processes) in ""''^^CaH-^''**Pb 
have shown different behaviours in the two systems. In par- 
ticular, at sub-barrier energies, reactions induced by '^''Ca were 
found to produce larger capture cross-sections [42]. This in- 
crease is consistent with our interpretation which is that this is 
a result of positive Q-value transfer reactions associated with 
N/Z equilibration in the ^°Ca reactions. 

To conclude, experimental MAD for reactions with small 
isospin asymmetry show that magic numbers in the entrance 
channel reduce quasi-fission and are thus expected to increase 
the probability for fusion, while non magic systems show more 
quasi-fission. With a large initial isospin asymmetry, a rapid 
N/Z equilibration occurs in the early stage of the reaction, mod- 
ifying the identities of the collision partners. This is the case for 
'**'CaH-^''^Pb, which, as far as the competition between fusion 
and quasi-fission is concerned, behaves more like a non magic 
system, i.e., with increased quasi-fission. Reactions with the 
neutron-rich ^^Ca on heavy targets usually have small isospin 
asymmetry, and thus are more favourable to fusion than reac- 
tions with ^°Ca, as well as leading to more neutron-rich com- 
pound nuclei having a higher probability of surviving fusion- 
fission. The importance of isospin asymmetry in the entrance 
channel should be considered in planning fusion experiments 
with exotic beams to form and study new isotopes of existing 
elements, as well as new super-heavy elements. 
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